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TARGET  PARAMETER  ESTIMATION  IN  THE  NEAR-FIELD  WITH  TWO  SENSORS 
HEMCHANDRA  M.  SHERTUKDE  and  YAAKOV  BAR-SHALOM 
Electrical  Engineering  and  Computer  Science  Department, 

University  of  Connecticut,  U-157, 

Storrs.  CT  06268 
ABSTRACT 

A  performance  prediction  procedure  is  developed  for  the  evaluation 
of  a  passive  tracking  technique  for  the  localization  of  targets  in  the 
near-field  of  the  sensors.  The  target  moves  in  a  three-dimensional 
space  with  a  constant  velocity  at  a  constant  depth.  The  parameter 
vector  is  thus  f  i  ve  -  di  me  ns  i  o  nal .  Target  localization  parameter 
identif  iabi  li  ty  is  established  with  the  aid  of  the  Fisher  Information 
Matrix  (FIM).  The  FIM  is  evaluated  for  various  combinations  of  the 
following  sequences  of  measurements:  time  difference  of  arrival  (TDOA) 
between  the  two  sensors,  depression  angles  measured  at  the  two  sensors, 
and  frequency  measurements  obtained  separately  with  the  two  sensors. 

The  FIM  is  used  to  determine  bounds  on  localization  performance, 
and  the  corresponding  uncertainty  ellipses  associated  with  the  target 
position  are  evaluated  for  various  tracking  scenarios. 

The  theoretical  results  are  corroborated  by  applying  a  maximum 
likelihood  estimation  algorithm  to  simulated  data  and  observing  the 

results  for  a  series  of  Monte  Carlo  simulations.  I  Accession  For 
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I.  INTRODUCTION 


The  problem  of  passively  tracking  the  position  of  a  target  in  two 
dimensions  has  been  extensively  discussed  in  the  literature  lll-{71.  A 
performance  prediction  procedure  is  developed  and  applied  to  the 
evaluation  of  a  passive  three-dimensional  tracking  problem,  where 
the  target  moves  at  a  constant  depth  with  constant  velocity.  The  position 
and  velocity  of  the  target  can  be  estimated  from  passive  measurements 
obtained  with  as  few  as  two  sensors  which  are  omnidirectional  in  the 
horizontal  plane  with  possible  added  capability  to  measure  depression 
angles,  in  which  case  they  are  vertical  arrays. 

The  Fisher  Information  Matrix  (FIM)  [10)  is  evaluated  for  various 
combinations  of  the  following  sequences  of  measurements  :  time  difference 
of  arrival  (TDOA)  between  the  two  sensors,  depression  angles  measured  at 
the  two  sensors,  frequency  measurements  obtained  separately  with  two 
sensors.  The  following  issues  are  addressed: 

1)  Identif iability  -  under  what  conditions  (if  any)  can  the  target 
localization  parameters  be  estimated. 

2)  Assessment  of  the  quality  of  the  parameter  estimates  via  evaluation 

of  the  Cramer-Rao  Lower  Bound  (CRLB)  of  the  estimation  error  covariance 
(91,  (101. 

3)  Development  of  a  suitable  numerical  algorithm  that  yields  the 


Maximum  Likelihood  Estimate  CMLE)  of  the  parameters  by  maximizing 
the  likelihood  function  [10]. 


4)  Comparison  of  the  performance  obtained  with  simulated  data  to  the 
theoretically  predicted  CRLB  for  validation  of  this  approach. 

The  results  demonstrate  that  reliable  estimates  of  the  target  position 
in  three  dimensions  and  its  velocity  can  be  obtained  with  measurements 
such  as  TDOA  between  the  two  sensors  and  two  depression  angle 
measurements.  Addition  of  frequency  measurements  enhances  the  accuracy  of 
the  estimates.  Even  with  a  very  limited  number  of  measurements,  e.g.,  six 
to  eight  TDOA  measurements  combined  with  the  depression  angle 
measurements,  the  target  can  be  localized  when  in  the  vicinity  of  the 
sensors,  i.e,  in  a  region  where  the  observed  measurements  change 
noticeably  with  time.  In  addition,  the  results  indicate  that  the 
near-field  estimates  obtained  from  combined  TDOA,  depression  angle  and 
frequency  measurements  are  comparable  with  bearing  measurements,  which 
require  additional  complexities.  The  table  below  summarizes  the  results 
for  various  sensor  configurations. 


Conf. 

TDOA 

Freq. 

Dep. 

Ang. 

Localization  performance 

A 

X 

- 

X 

Identifiable 

B 

X 

- 

- 

Practically  unidentifiable 

C 

X 

X 

X 

Identifiable  (good  performance] 

D 

- 

X 

- 

Identifiable  (marginal  performance) 

E 

X 

X 

- 

Identifiable  (better  than  D) 

F 

- 

X 

X 

Identifiable 

G 

- 

- 

X 

Practically  unidentifiable 

Table  1.  Summary  of  results 


Section  II  discusses  the  formulation  of  the  target 
localization  problem  in  three  dimensions,  and  the  i de n t  i  f  i  a b i  1  i  t y 
condition,  stated  in  terms  of  the  Fisher  Information  Matrix. 

In  addition,  other  measures  of  performance  relating  to  eigenvalue 
dispersion  and  target  position  uncertainty  are  also  discussed. 

Section  III  deals  with  the  CRLB  for  target  localization  which  is 
numerically  evaluated  for  various  target  trajectories  and  sequences  of 
measurements. 

In  Section  IV  the  target  localization  performance  is  evaluated 
experimentally  using  simulated  data  in  conjunction  with  a  gradient-based 
Maximum  Likelihood  algorithm  and  the  errors  are  compared  to  the  CRLB.  The 
results  compare  well  with  the  theoretical  predictions.  Examination  of  the 
performance  for  various  test  cases  provides  important  insight  into  the 
localization  accuracy  as  it  depends  on  the  target  trajectory  relative  to 
the  sensor  positions  and  the  types  of  measurements  used. 

Section  V  presents  a  summary  and  conclusions.  The  approach  is 
sufficiently  general  to  be  applied  to  air  acoustics,  passive  radar  or 
electromagnetic  direction  finding. 


11.  PARAMETER  I  DENT  IF  I  ABILITY  AND  PERFORMANCE  PREDICTIONS. 


The  target  localization  information  contained  in  sequences  of 
measurements  is  examined  with  the  aid  of  the  Fisher  Information  Matrix. 

The  FIM  will  be  used  to  determine  the  CRLB  on  the  error  covariance  matrix 
of  the  target  localization  parameter  vector  6. 

The  analysis  is  performed  for  a  three-dimensional  case  where  the 
target  trajectory  and  the  sensors  are  not  in  the  same  plane.  In 
configuration  A  from  Table  1  the  passive  acoustic  sensors  (vertical 
arrays,  omnidirectional  in  the  horizontal  plane)  are  placed  on  the 
surface  of  the  ocean  and  the  target  moves  at  a  constant  depth. 

The  analysis  is  done  for  various  combinations  of  the  following  sequences 
of  measurements  : 

1)  TDOA  between  the  two  sensors  :  T^(0,t^). 

2)  Depression  angle  from  sensor  i  :  AjO.t^,),  i-1,2 

3)  Frequency  (Doppler-shifted)  at  sensor  j  :  fj(0,t^),  j“l,2 

II. 1  CRLB  for  TDOA  and  Depression  Angle  Measurements. 

Consider  the  two  sensors  represented  in  Fig.l  with  known  locations 
(  Xj  ,  y^  ,  0  )  ,  i-1,2.  The  position  of  the  target  at  time  t^  is  given  by: 

(11.1.1) 

(11.1.2) 


x(tk)-x„*v,t^ 

y(t^)-y^-Vyt^ 

z(t^)-z„ 


(II.1.3) 


The  parameter  vector  9,  describing  the  initial  position  and  velocity,  is: 

Vo*  ''x>  V  (II. 1.4) 

where  v, ,  Vy,  denote  the  x  and  y  components  of  the  constant  velocity 
vector  of  the  target,  and  '  denotes  vector  or  matrix  transpose.  The  target 
moves  at  a  constant  depth  z^. 

Assuming  the  TDOA  estimates  (measurements)  Tj(6,t^)  are 
proportional  to  the  corresponding  range  difference  rj(0,t^)  one  has 

'^d=''d(0.tkVc  (II.1.5) 

where, 

rd(0.tk)-r,(0,t^)-r2(9,^)Sh,(e,^) 

1 

I 

-{ CIU.6) 

and  c  denotes  the  speed  of  sound  propagation  In  the  medium.  The 
corresponding  range  difference  measurement  Zj(t^)  made  with  a  sensor- 
pair  is  given  by 

z,(t^)-h,(0,t^)-Wj^(^),  k-l,.,n  (11.1.7) 

where  w  denotes  the  measurement  noise. 

The  depression  angles  are 

h2(e,tk)SA,(0,tk) 


[(x,-x(t^))2*(yj-y(t^))2]2 


(II.1.8) 


and. 


h3(0,t^)^A2(e.^) 


[(X2-x(y)2*(y^-y(t^)}2]2  j 

The  correepondine  measurements  made  with  the  sensors  are 

Z2(tk}-h2(0.tk)+w,^(tk).  k-l...n  (II.l.lO) 

and, 

Z3Ct^)-h3(0.t^)+w2^(t,^),  k-l..,n  (II.l.ll) 

where  w^^  ,  w^^  denote  the  angle  measurement  noises. 

Assuming  that  the  measurement  noises  are  Independent,  identically 
distributed,  normal  with  zero  mean  and  variances  ®z,^  ■ 

i-1,2,3,  the  individual  likelihood  functions  of  the  parameter  0  (the 
Joint  pdf  of  the  observations  Z,  =  [z,  (t,^),  k»l,..,n],  1-1, 2, 3, 

conditioned  on  0)  can  be  written  as: 

A2j(0)-p(2,|0) 

■  ft  p(Z,(U)|0) 
k-1 

=  «.exp  ^  §  [  z,(tk)-h,(0,t^)  ] ^  j.  (II.1.12) 

where  i  =  l,2,3  and  k=l,  .  .,n,  and  a,  denotes  the  normalization 

constants. 

A 

The  covariance  of  any  unbiased  estimate  0(Z)  of  the  (non- 
random)  parameter  0  cannot  be  smaller  than  the  CRLB  namely 

P  -  E  [(0-0)(0-0)-|0]  i  [Jzf®)]  *  (II.1.13) 


where 


(II.1.14D 


is  the  FIM  corresponding  to  the  combined  informational  content  of  the 
three  independent  sets  of  measurements.  The  Hessian,  symbolically  written 
above  as  the  second  derivative,  is  evaluated  at  the  true  value  of  the 
parameter  0.  The  additivity  follows  from  the  assumption  of  Independence 
of  the  measurement  noises. 

It  can  be  seen  that  the  i  de  n  t  i  f  i  ab  i  li  t  y  condition  is  equivalent  to 
requiring  the  FIM  to  be  nonsingular.  Otherwise  the  lower  bound  is  not 
finite  and  the  parameters  are  not  identifiable. 

Using  the  gradient  notation  (a  column  vector)  where 


V  _d_  _d_  ^“l' 

0  |_  dx„’  dy^'  dv,’  dVy'  dz„  J 

one  can  evaluate  (II. 1.14)  as  follows; 

Jz(0)  -  §  -  EfVeVe'lnA,  (0)] 

1-1  ‘ 

i-1  *1  'k-1 

3  (  n 

E  Z  (V,h,(0,L  )](Veh,(0,t^))- 
1-1  '  ^k-1 

-(7,V,-h,(0,t^)]{z,(t,^)-h,(0,t^)]J 


(11.1.15) 


(II.1.16) 


Denoting  the  gradient  and  Hessian  of  h  as  hg  and  hg^ , 
respectively,  i.e.. 


(II. 1.17) 


V9V9-h,(0.t^)=h,ee(tk)  (II.1.18) 

for  i  -  1 , 2 , 3  yields 

Jz(0>|;-izi(®i 

i-1 

3  ,  n  V 

-5;o^-2  E{2h,ettJhie'^tk)-h,ee(tJ[z,(t^)-h,(04)l}  (11.1.19) 

1«1  ‘  k=l 

Since  expressions  (II. 1.17)  and  (II. 1.18)  are  d  te  r  m  i  n  i  s  t  i  c  ,  the  only 
stochastic  term  in  (II. 1.19)  is: 

z,(tk)-h,(04)=  w^^(tk)  (II.1.20) 

which  was  assumed  to  be  zero-mean.  Thus  (II. 1.19)  becomes 

Jz(0>i°z?i»^l9ttv.5hi9(t^)  (11.1.21) 

I»1  ‘k=l 

where  it  is  assumed  that  the  set  of  TDOA  and  depression  angle  measurements 
are  obtained  with  a  single  sensor-pair  at  n  discrete  instants  of  time. 

A  necessary  condition  for  the  invertibility  of  (11.1.21)  is  that  n2dim(0)-m, 
in  this  case  m-5.  The  analytical  expression  for  (II. 1.21)  can  be  obtained 

from  (II. 1.6),  (II. 1.8),  and  (II. 1.9)  but  at  the  same  time  it  is  too 
complicated  to  obtain  explicit  expression  for  its  invertibility.  As 
discussed  in  [8]  the  full  rank  condition  of  the  set  of  gradients  in 
(II. 1.21)  Is  sufficient  for  Invertibility  and  thus  local  1  de  n  1 1  f  i  ab  1 1 1 1  y . 
Since  our  problem  pertains  to  local  identifiability  one  can  work  on  this 
condition  to  obtain  the  results. 

The  special  structure  of  (II. 1.21)  as  a  sum  of  dyads  Implies  that  it 
will  be  invertible  iff  the  set  of  vectors  h,9(t^),  k-l,.,n,  i-1, 2, 3, 


span  the  Euclidean  space  of  dimension  m-dim(6). 


II. 2  The  Identif iability  Condition 

To  determine  target  localization  parameter  identifiability  from  TDOA 
and  the  two  depression  angle  measurements,  the  gradients  hjgCt,^),  -1,2,3  are 
evaluated  using  (11.1.6),  (II. 1.8)  and  (II. 1.9). 

We  denote  the  gradient  vectors  as  follows 


C  ®ik*  ^ik«  ®lk^»  ^ik^*  ‘'ik  3  •  i-1,2,3  (II.2.1) 

where  the  quantities  in  the  last  bracket  can  be  obtained  explicitly  as 
done  in  Appendix  A.  When  only  TDOA  measurements  are  considered  the 
spanning  condition  is  equivalent  to  the  determinant  being  non-zero, 
where 


ID. 


®12 

*13 

®l-4 

®15 

bit 

b.3 

b.5 

®ll^l 

®i2^2 

®i3^ 

®M^-I 

®J5^5 

bi,tt 

^.2^2 

b,3t3 

bHt„ 

b,5t5 

‘'ii 

^i2 

Ci3 

X  0,  i-1 


(II.2.2) 


Similarly  I D2 1  and  I  D3 1  can  be  written  corresponding  to  depression  angle 
measurements  by  substituting  i-2,3  in  Eq.  (II. 2. 2),  respectively.  However, 
if  all  the  sets  of  measurements  are  also  considered  the  determinant  of 
(II. 1.21)  should  be  considered.  Since  the  determinant  of  the  sum  of 
matrices  is  very  difficult  to  evaluate  the  identifiability  will  be 


discussed  in  the  sequei  for  one  type  of  measurement  at  a  time.  The 
evaluation  of  the  determinants  is  relegated  to  Appendix  B. 

The  conclusions  arc  summarized  next.  When  the  target  is  at  endfire, 
i.e.,  target  bearing  angle  is  zero  and  target  depression  angles  are  also 
zero,  then  the  gradients  are  zero,  which  precludes  i de n t  i  f  i  a b i  1  i  t y . 
Furthermore  in  the  far-field  B,(0,t^)SB2(0.t^)  and  Aj(0,t^)=A2(0.t^)  which  also  precludes 
identifiability.  When  B,  C0,t,j)'*-B2(0.t^)-n,  and  A,  (0,t^)-A2(0,t^)  i.e.,  target 
motion  in  the  perpendicular  bisector  plane  of  the  segment  joining  the 
sensors  one  has  bj^-0,  Cy^-O  and  the  target  cannot  be 
practically  Identified  as  will  be  seen  later.  From  the  above  the 
near-field  condition  is  necessary  for  identifiability.  The  condition  for 
|D,i><0  has  been  discussed  in  [7]  for  a  2-d  case.  To  obtain  a  similar 
condition  for  a  3-d  the  case  is  tedious  and  has  been  relegated  to 
Appendix  B.  Thus,  with  the  5-dimensional  parameter  vector  for  a  target 
moving  in  a  3 -dimensional  space  one  has  practical  Identifiability  from 
TDOA  and  depression  angle  measurements  unless; 

1)  Target  at  endfire  or  far-field. 

2)  Target  motion  in  the  perpendicular  bisector  plane  of  the  segment 
Joining  the  two  sensors. 

II. 3  CRLB  for  the  Combination  of  Complementary  Measurements. 

The  contribution  of  additional  measurements  made  with  a  pair  of  sensors, 
such  as  single-sensor  target  frequency  estimates  (measurements],  which  are 


denoted  with  i-4.5  as 


fj(e.tk)feh,(0,tk)-(l-r,_3(e.t^)/clfo  i-4,5  j-1,2  (II.3.1) 

where  r,_3(0,t^)  denotes  the  time  derivative  of  the  range  measured 
by  the  respective  sensor,  c  denotes  the  speed  of  sound  propagation  in  the 
medium  and  fo  denotes  the  Doppler-free  target  frequency.  The 
corresponding  measurements  are 

z,(tk)-h,(0.t^)^w^^(t^)  i-1.2  (II.3.2) 

where  w,  is  the  measurement  noise,  which  is  white,  with  zero  mean 

and  a  variance  a,  ^  . 

The  total  FIM  Jz(0),  corresponding  to  all  the  sets  of  measurements 
obtained  with  a  pair  of  sensors,  is  given  by 

Jz(9)-  2  h,^(t^)h,e'(tk)  CII.3.3) 

i-1  '  k-1 

where  hi^ft^)  denotes  the  gradient  of  h^fO.t^)  as  defined  earlier  in  Eq.  (II. 1.17). 
//.4  Additional  Performance  Measures. 

The  inverse  of  FIM  is  the  CRLB  for  the  error  coavriance  matrix  P 
associated  with  the  target  localization  parameter  estimates  i.e., 

P2J"‘(0)  (II.4.1) 

where  J“‘(0)  denotes  the  FIM  for  the  combination  of  multiple  measurement 

sets.  The  determinant  of  J"’(0)  provides  a  measure  of  the  lower  bound 
of  the  hypervolume  of  the  multiple  dimensional  uncertainty  ellipsoid  [10], 
associated  with  target  localization  parameter  estimate.  The  shape  of  the 


uncertainty  ellipsoid  can  be  obtained  from  the  solution  of  the  eigensystem 


involving  The  eigensystem  provides  the  following  interpretations: 

i)  The  lengths  of  the  axes  and  their  orientation  in  the  reference  coordinate 
system  are  given  by  the  corresponding  eigenvalues  and  eigenvectors,  i.e., 

J'‘(0)E.  -  KjEj  (II.4.2) 


where  Ej  and  are  the  eigenvectors  and  the  corresponding  eigenvalues,  respectively, 
ii)  The  dispersion  of  the  eigenvalues,  ^n>ax  ^  ^min  •  also 
provides  a  measure  of  parameter  i  den  tif  iab  i  1  i  ty  [7], 

The  uncertainty  ellipse  for  the  target's  horizontal  position  at 

time  t  can  be  easily  derived  from  J”'{0)  as  follows.  The  predicted 
position  for  time  t  is 

[Jffi]-[J  f  i  ?  8]«  ‘“•‘>■3) 

One  can  obtain  the  corresponding  position  error  covariance  matrix 
corresponding  to  time  t,  Cp(0.t)  as: 

[if]  CII.4.4) 

The  shape  and  area  of  this  position  uncertainy  ellipse  can  be  obtained 
from  the  appropriate  eigensystem  solution. 


III.  PERFORMANCE  PREDICTIONS 


The  performance  measures  discussed  in  Section  II,  i.e., 
identi f  iab il i  ty  .  matrix  or  data  ill -conditioni  ng  ,  the  area  and  elongation 
of  the  ellipse,  are  determined  for  various  situations  in  this  section. 

Five  scenarios  of  target  trajectory  and  sensor  configurations  were  chosen 
to  determine  the  feasibility  of  target  localization  using  measurements 
from  a  pair  of  sensors.  Using  sequences  of  TDOA,  depression  angles  and/or 
frequency  measurements,  questions  of  identif  i  abi  lity ,  matrix 
ill-conditioning,  and  target  parameter  uncertainty  will  be  addressed  with 
the  aid  of  FIM. 

Case  I 

The  first  case  consists  of  a  target  trajectory  making  an  angle  of 

45°  with  the  array  axis  as  shown  in  Fig. 2a.  The  target  passes  between 
the  sensors  spaced  2  km.  The  target  location  in  this  case  is  generally 
identifiable,  from  TDOA  and  depression  angle  measurements  (sensor 
configuration  A  from  Table  1).  The  FIM  was  determined  for  a  constant  speed 
of  10  kn  using  a  sequence  of  eight  TDOA  and  depression  angle 
measurements  corresponding  to  the  eight  equispaced  target  positions 
indicated  in  Fig. 2a.  The  target  is  assumed  to  be  moving  at  a  depth  of 
0.2  km.  The  FIM  was  evaluated  and  found  to  be  non-singular  with  a 
conditioning  measure,  log^jj(k„„/k„,„)«7.7  8.  The  standard  deviations  of  TDOA  and 

depression  angles  are  0.1s  and  6°  respectively.  The  area  of  the 
corresponding  two-sigma  horizontal  position  uncertainty  ellipse,  shown  in 

Fig. 2a  corresponding  to  the  largest  ellipse  is  0.9  k  m^.  Under  the 
Gaussian  assumption,  the  positional  estimate  will  fall  within  this  ellipse 


with  86.5  %  probability.  The  depth  standard  deviation  was  0.04  km.  The 
addition  of  frequency  measurements  to  the  above  measurements 
(configuration  C)  further  improves  the  performance.  The  zero-Doppler 
frequency  is  100  Hz.  The  standard  deviation  of  frequency  measurements  was 
O.lHz.  The  corresponding  conditioning  measure  was  7.55  and  the  largest 
area  of  the  corresponding  two-sigma  position  uncertainty  ellipse  was  0.615 
km^  as  shown  in  Fig.  2b  and  the  depth  standard  deviation  was  0.0358  km. 
Using  only  TDOA  and  depression  angle  measurements  with  the  latt  ,  more 

accurate  to  a  =  3°  resulted  in  a  considerable  improvement  in  ‘he 
performance  measure.  The  corresponding  values  for  the  conditioning  measure 
and  for  the  area  of  the  largest  two-sigma  uncertainty  ellipse  are  6.96  and 
0.319  km^  respectively.  The  results  are  shown  in  Fig. 2c  and  the  depth 
standard  deviation  was  0.0348  km. 

It  is  interesting  to  note  that  results  from  Fig. 2a  when  compared  to 
those  from  (7]  show  that  for  these  numbers  the  3-d  problem  with  TDOA 
and  depression  angles  yields  better  accuracy  than  the  2-d  problem  with 
TDOA  measurements  only. 

Case  II 

The  second  case  examined  consists  of  a  target  trajectory  that  is 
perpendicular  to  the  array  axis  and  intersecting  the  axis  at  a  point 
outside  the  sensors  as  shown  in  Fig. 3a.  Configuration  A  namely,  TDOA  and 
the  two  depression  angle  measurements  was  considered.  The  FIM  was  very 
ill-conditioned  and  the  uncertainty  areas  are  too  large.  The  addition  of 
frequency  measurements  (configuration  C)  significantly  improves  target 
localization.  The  FIM  was  first  obtained  for  the  TDOA  and  the  two 
depression  angle  measurements  (configuration  A),  each  set 


consisting  of  eight  samples.  For  a  target  speed  of  10  kn  and  the  sensors 
spaced  2  km  apart,  the  matrix  conditioning  measure  is  7.18  and  the  area 
of  the  two-sigma  position  uncertainty  ellipse  shown  in  Fig. 3a  is  1.6  km^ 
and  the  depth  standard  deviation  was  0.0634  km.  The  addition  of  two  sets 
of  frequency  measurements  and  all  the  other  quantities  remaining  the  same 
(configuration  C),  gave  a  matrix  conditioning  measure  of  6.48,  and  the 
area  of  the  two-sigma  position  uncertainty  ellipse  as  0.496km' 
respectively,  as  shown  in  Fig. 3b  and  the  depth  standard  deviation  was 
0.05698  km. 

Case  III 

The  third  case  examined  consists  of  a  target  trajectory  that  is 
parallel  to  the  array  axis.  Evaluation  of  the  FIM  based  on  the  TDOA 
and  the  depression  angle  measurements  (configuration  A)  indicated  still 
poor  localization  as  in  Case  II.  This  is  illustrated  in  Fig. 4a. 
Localizability  can  be  improved  by  the  addition  of  the  frequency 
measurements  (configuration  C).  The  FIM  was  obtained  for  combined 
sets  of  TDOA  and  depression  angle  measurements  consisting  of  eight 
samples  each  and  sets  of  TDOA,  depression  angles  and  frequency 
measurements,  consisting  of  eight  samples  each.  The  matrix  conditioning 
measure  and  area  of  the  two-sigma  position  uncertainty  ellipse  as  shown  in 
Fig. 4a.  and  Fig. 4b  are  7.43,  and  1.357  km^  (configuration  A) 
respectively,  and  5.98  and  0.438  km^  (configuration  C)  respectively, 
and  the  depth  standard  deviations  were  0.099  km  and  0.089  km. 


The  fourth  case  consists  of  a  target  trajectory  perpendicular  to  the 
array  axis  passing  through  the  individual  sensors  and  equidistant  to  the 
two  sensors  spaced  at  2  km  i.e.,  in  the  perpendicular  bisector  plane. 

The  FIM  was  evaluated  on  the  basis  of  TDOA  and  depression  angle  estimates 
(configuration  A].  Regardless  of  target  initial  position  along  this 
trajectory  or  target  speed,  a  constant  TDOA  sequence  of  zero  results.  The 
target  localization  parameters  are  practically  unidentifiable  if  based  on 
TDOA  and  depression  angle  measurements  alone,  as  shown  in  Fig. 5a.  To 
localize  a  target  with  this  trajectory,  it  is  necessary  to  obtain 
additional  measurements.  Frequency  measurements  obtained  with  each  of  the 
sensors  were  used  in  conjunction  with  the  TDOA  and  depression  angle 
measurements  (configuration  C).  The  FIM  was  evaluated  for  the  combined 
five  seta  of  measurements.  The  matrix  conditioning  measure  is  6.73, 
which  although  relatively  large,  produces  resonable  results.  The 
corresponding  two-sigma  position  uncertainty  ellipse  illustrated  in  Fig. 5b 
has  an  area  of  0.205km^  and  the  depth  standard  deviation  was  0.0367  km. 
The  FIM  was  also  evaluated  for  depression  angle  measurements  only 
(configuration  G)  and  the  target  was  found  to  be  practically 
unidentifiable.  The  standard  deviation  of  the  depression  angle 

measurements  was  6°.  The  conditioning  measure  was  6.73  and  the  area  of 
the  two-sigma  uncertainty  ellipse  was  17.41  km^  as  shown  in  Fig. 5c  with 
the  depth  standard  deviation  of  0.0583  km.  Finally  the  FIM  was  evaluated 
for  frequency  measurements  only  obtained  separately  with  the  two  sensors 
(configuration  D).  The  standard  deviation  of  the  frequency  measurements 
was  O.lHz.  The  resulting  conditioning  measure  was  7.45  and  the  area  of  the 
two-sigma  uncertainty  ellipse  was  2.59  km^  as  shown  in  Fig.Sd,  while  the 


depth  standard  deviation  was  4.46  km.  The  horizontal  localization 
uncertainty  is  not  too  large  but  this  measurement  carries  practically  no 
information  about  depth. 

Case  V 

In  this  case  the  target  trajectory  was  as  in  Case  I  but  the  types  of 
measurements  are  different.  The  FIM  was  evaluated  using  only  the 
target  frequency  measurements  obtained  separately  with  two  sensors 
(configuration  D).  The  results  revealed  that  the  target  is  generally 
identifiable  as  shown  in  Fig. 6a.  Addition  of  depression  angle  measurements 
(configuration  F)  improved  performance  as  shown  in  Fig. 6b.  The 
conditioning  measure  and  the  area  of  the  largest  two-sigma  uncertainty 
ellipse  are  7.097,  and  0.805  km^,  respectively,  for  (configuration  D) 

and  the  depth  standard  deviation  was  2.16  km;  and  6.74  and  0.325  km^, 
respectively,  (configuration  F)  as  in  Fig. 6b.  and  the  depth  standard 
deviation  was  0.0384  km  respectively.  Finally  the  FIM  was  evaluated  by 
assuming  only  TOOA  and  frequency  measurements  (configuration  E).  The 
target  is  identifiable  and  the  corresponding  matrix  conditioning  measure 
is  7.467  and  the  area  of  the  two-sigma  position  uncertainty  ellipse  as  in 
Fig. 6c.  is  0.203  km^  and  the  depth  standard  deviation  was  1.211  km. 

Note  that  in  all  the  above  cases  that  have  been  evaluated  the  target's 
direction  of  travel  along  the  constant  velocity  trajectory  does  not  affect 
the  results  of  the  analysis. 


IV.  VALIDATION  OF  THE  PERFORMANCE  PREDICTION  TECHNIQUE. 

Experimental  performance  evaluation  was  done  with  the  implementation 
of  a  maximum  likelihood  parameter  estimation  algorithm  by  using 
simulated  measurement  data.  The  algorithm  uses  the  Conjugate  Gradient 
Search  method  to  minimize  the  least  squares  criterion  equivalent  to  the 
maximum  likelihood.  In  this  algorithm  the  error  measure  is  the  difference 
between  the  noisy  measurement  data  and  measurement  predictions  based  on 
the  estimated  localization  parameter  values.  Simulated  data  consisting  of 
TDOA,  depression  angle  and  frequency  measurements  were  formed  for  one 
selected  scenario.  The  measurement  data  were  corrupted  by  additive  white 
Gaussian  noise  of  a  known  variance.  For  one  scenario  as  in  Case  I,  50 
Monte  Carlo  tracking  simulations  were  made  to  form  a  statistical  base  for 
examination  of  the  localization  error.  The  resulting  localization 
estimates  were  then  plotted,  forming  scatter  diagrams  indicative  of  the 
algorithm's  performance.  The  corresponding  two-sigma  error  ellipse  as 
predicted  by  the  FIM's  were  superimposed  over  the  scatter  plots,  allowing 
for  a  direct  comparison  between  the  experimental  and  the  theoretically 
predicted  performance.  The  results  are  shown  in  Fig. 8.  For  the  50  tests 
performed,  good  agreement  exists  between  the  experimental  results 
and  the  theoretical  uncertainty  ellipses  shown.  A  sequence  of 
uncertainty  ellipses  and  the  corresponding  scatter  diagrams  are  shjwn  in 
the  figure  for  the  duration  of  the  target  path  from  which  measurements 


were  obtained. 


V.  SUMMARY  AND  CONCLUSIONS 


Following  the  earlier  work  [7],  where  the  analysis  was  done  for  a 
two-dimensional  tracking  problem  with  two  omnidirectional  sensors, 
in  this  work  the  analysis  has  been  extended  to  a  three-dimensional 
tracking  problem  with  two  sensors.  An  analysis  procedure  based 
on  the  information  contained  in  various  sets  of  measurements,  viz.  the 
TDOA,  the  two  depression  angles,  and  the  target  frequencies  measured 
separately  with  the  two  sensors  has  been  developed  and  tested  for 
evaluation  of  target  localization  with  a  pair  of  sensors,  when  the  target 
moves  in  the  vicinity  of  the  sensors.  This  tool  was  initially  used  to 
determine  if  the  informational  content  of  a  set  of  measurements  is 


sufficient  to  identify  the  localization  parameters.  The  analysis  also 
yields  important  performance  measures  such  as  algorithm  convergence 
indicators,  position  and  velocity  parameter  uncertainty  ellipses  for  a 
three-dimensional  passive  tracking  problem.  Examination  of  the  performance 
for  various  test  cases  provided  insight  relating  localization  accuracy  to 
target  trajectory  relative  to  sensor  positions,  measurement  types  and 
measurement  quality.  The  results  demonstrate  that  reliable  localization 
based  on  TDOA  and  depression  angle  measurements  occurs  primarily 
for  targets  passing  between  the  sensors.  The  TDOA  measurements  in 
conjunction  with  the  depression  angle  measurements  were  observed  to  yield 
a  somewhat  elongated  position  uncertainty  ellipse.  The  simultaneous  use  of 
all  types  of  measurements,  i.e.,  frequency,  TDOA,  and  depression  angle 


produced  an  improved  uncertainty  ellipse  in  terms  of  eigenvalue 
dispersion. 


A  more  interesting  finding  is  that  a  target  in  motion  in  S-dimensions 
can  be  localized  by  using  only  the  target  frequency  measurements  obtained 
separately  with  two  sensors. 

The  results  of  the  computer  simulations  are  in  good  agreement 
with  the  theoretical  performance  predictions  obtained  with  the  CRLB. 

The  following  important  results  can  be  stated  : 

1)  A  target  in  motion  in  three  dimensions  can  be  localized  by  evaluating 
the  TDOA  and  depression  angle  measurements  (configuration  A)  obtained  with 
two  sensors,  unless  the  target  is  at  endfire  or  far-field.  Motion  in  the 
perpendicular  bisector  plane  is  practically  not  identifiable. 

2)  A  target  in  motion  in  3-dimensions  cannot  be  practically  localized  by 
evaluating  only  the  TDOA  measurements  (configuration  B)  obtained  with  two 
sensors.  However,  it  can  be  localized  by  evaluating  only  the  target 
frequency  measurements  obtained  separately  with  two  sensors 
(configuration  D).  Any  additional  measurements  like  the  TDOA  or  the 
depression  angle  measurements  improve  performance. 

The  analysis  is  general  and  can  be  applied  to  air  acoustics,  passive 
radar,  and  electromagnetic  direction  finding.  With  such  an  analysis  tool 
one  can  predict  performance,  i.e.,  the  ability  to  track,  and  separate 
targets,  for  various  types  and  qualities  of  measurements.  This  analysis 
Indicates  the  potential  of  a  field  of  sensors  in  tracking  targets. 
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APPENDIX  A 


Evaluation  of  Gradients  - 

To  determine  target  localization  identif i abi li ty  from  TDOA  and  the 
depression  angle  estimates,  the  gradients  of  hj  ( 0 ,  t^^ )  ,  i-1,2,3 

evaluated  using  equations  (II. 1.6),  (II. 1.8)  and  (II. 1.9)  in  Section 
as  folio  WS: 

d  .  ^  x,-x(t^)  X2-x(t^) 

/x(t^)-xa  /p,ce,t,)|  /X(tk)-X2^ 

°  l^Q.tkU  ■  IpT^I 

*C0SB^C03A^Jj  -  C03B2kC0SA2t 
“  ®tk 

where 


P,C0.tk)= 


[ 


(K^-x(tfc))2^(y,-y(tk))2 


Y 


and  , 


i 

P2Ce.t^)=  [  (X2-x(tk))2*(y^-y(t^))2]2 


Ik 


7k 


two 

are 

II. 


(A.1) 

(A.  2) 

(A.3) 


and  B 


and  B 


denote  the  bearing  angles  at  time 


measured  with  respect  to  a  plane  passing  through  the  two  sensors 


and  the  target.  Similarly 


^h,(0,t^)  -  sinB^^cosAy^  -  sinB2^cosA2k 


-  sinA^  -  sinA^j 


-  c 


Ik 


y 


Similarly  gradients  of  hjfB.tjj)  and  bgfS.tjj)  are 


-r - ^ - ii 

\  [(xrxdiJl'-lyryOi')']  ) 


f  Zq  ]  f(x,-xCtfc))> 
’  lr,(0.t,)2/  I  p^te.t,)  / 


sinA, 


.-.■„^(_C03B^^) 


4:*--  r - ^ - r  1 

^  [<x,-x(tj))^*(y,-yCt,))']^  J 


sinA,. 

=— rT^C-sinB,,^) 


(A.  4 


(A.5 

(A.6 


(A.7 


(A.8 


dv 


dv  "  ’^2k^k 

y 


^J^hjCe,.^)  -  {  °  2}  (p'^ie.'v)} 


-rj-(-cosB2^) 


■  {r,u°t,A  {p/e.v} 


.rjte.tfc) 


=— r^C-sinBjk) 


cosAjjj 

is 


dVjj*^3^®''t^  "  ®3k^k 

^^hjCe.l,^)  - 

y 

If  the  frequency  measurements  are  also  considered  the  gradients  of 
i-4,5  are  evaluated  using  the  Eq.  (II. 3.1),  i.e. 

j^.ce.tk)  ■  -  -T['',(Xi-x(k))^*Vy(y^-y(t,^)(Xj-x(tk))]} 


8] 


'b.k.  i-4.5;  J-1,2 


^.CO.t,)  -  ^ - -3— i - 


i=4.5  j=1.2 


^hi(e.t^)  -  -c[ - f] 

[v,(xj-x(t^))*VyCyj-yCt^))  ](Xj-xCtt))t^i 


-dik.  i-4,5:  j-1,2 


J_h  re  t ,  - 

[v,(xj-x(^)j*VyCyj-y(t^))](yj-yCt^))t^i 


(A.19) 


(A.20) 


(A.2n 


.V/A 

v^>:c 


i-4.5i  j-1.2 


(A.22) 


Thus  the  gradient  vectors  h,g  (  0  ,  t,j )  ,  i-1,2,3,4,5  are  given  by: 


h,(0.tk)-Veh,(e.t^) 


C  ®lk' 


^Ik*  ®lk^'  ^Ik^t*  ‘'IkD 


(A.23) 
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APPENDIX  B 


Identif (ability  Conditions  for  the  3-d  case 
The  iden  ti  f  i  abi  ii  ty  condition  is  evaluated  in  two  steps: 

1)  To  prove  the  ide  nti  f  i  ab  i  1  i  ty  of  the  corresponding  system  only  for 
TDOA  estimates. 

2)  Apply  the  approach  in  step  (1)  to  the  general  problem  where  the 
TDOA  and  depresssion  angle  estimates  are  considered.  [In  addition 
the  complementary  measurement  sets  may  also  be  considered.) 

B.l  Step  !-  only  TDOA  measurements  considered. 

In  this  case  one  has  now  to  show  that  the  determinant  of  the  matrix 


D|  is  non-zero. 

Using 

the 

relation 

tfc  -  (k-1 )A 

• 

ID, 

1  reduces 

as 

follows: 

=t2 

ai3 

®15 

^12 

b,3 

b,5 

ID,|  - 

^11 

^12 

^t3 

^15 

(B.l. 

0 

®12^ 

a,32A 

0 

b,2A 

b,32A 

This  determinant 

can  be  developed 

along 

the  first  column 

from  which 

one 

can  deduce  that 

1  1  can  be 

zero  if 

1)  a,,  »  0  and,  (B.l. 2) 


2 )  b,, 


0  and, 


(B.1.2 


which  further  implies  that 


1) 

aik  ■  0 

and. 

(B.1.5) 

2) 

btt  -  0 

and. 

{B.1.6) 

3) 

o 

1 

O 

(B.1.7) 

{This  follows  from  the  fact  that  I  Dj  |  can  be  developed  along  any  column. 
Also,  since  the  minimum  number  of  measurements  are  considered.} 

«  0  =»  cosByjCosAyj  -  cosB2|jCOsA2k  •  0  (B.1.8) 

b^  -  0  =»  sir^^jjCOsAj^  -  sinB2,jC03A2k  =  0  (B.1.9) 

c^  -  0  sinAyj  -  sinAjij  -0  (B.1.10) 

The  above  three  conditions  imply  the  following: 
When  B^C0,tjj)*B2(0,tjj)  and  A,  (0,t,t)  =  A2(0,t^)  one  has  h,(0,t^)  *  0. 

Thus  a  measurement  does  not  provide  any  information,  when  the  target  is  at 
endfire.  Furthermore,  in  the  far-field  B,  (0,t,j)3^B2(0,t^) 
and  Af(0,t^)^A2(0.t,c)  and  all  the  gradients  are  practically 
zero,  which,  as  intuitively  expected,  precludes  identifiability.  In  other 
words  the  near-field  condition  is  necessary  for  identifiability.  When 
B^(0,t^)*B2(0,t^)  -  n,  and  A^(0,t^)  -  AjlO,!^)  i.e.,  along  the 

perpendicular  bisector  plane  of  the  line  joining  the  two  sensors,  one 
has  bf,^  0  and  Cy^  »  0  and  motion  along  this  perpendicular 

bisector  plane  also  precludes  identifiability. 

In  addition  for  determinant  of  0^  to  be  zero  the  following  have  to 

be  simultaneously  zero  i.e.. 


4)  ID^I 

-  0. 

and 

(B.1.11) 

5)  1D.J 

-  0. 

and 

(B.1.12) 

6)  1D.J 

•  0 

(B.1.13) 
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For  conditions  (B.1.11)  through  (B.1.13)  to  be  true  one  has  to  examine  the 

determinants  individually.  By  properly  partitioning  one  can  derive  the 

following: 


“  Isbell  I  I ^I>c22”^l>c2i^bcii  '^bc^l  °bcii  '  exists, 
where  it  can  be  shown  that 


|D»,„|.4co9^,cos(^^l^)sin(^ii2^)cos(^a^)^^^ 


-  4.cos»^co3(^'^^^^)sin(*'^^^^^^)co3(^*^^^^^)sin(^^^^^) 


(B.14) 


VJ.v.VJAV'V  ■■■ 


B-4 

where  ,C,  (0,t^)  and  C2(0,tjj)  are  angles  as  defined  in  the 

Fig.  B.l.  The  above  two  terms  are  identically  zero  iff: 

^  ^  C|2*^22*^ 

ii  )  C|2^C22*^ 

=>  Ca,*C2k-n  (B.1.15) 

i.e.,  motion  along  perpendicular  bisector  plane  of  the  line  Joining  the 
two  sensors. 

iii)  <t>^  -  J 

iv)  A,2-A22-n 

v)  A,3-'A23-n 

=*0k“2^  and  A„j^A2,j-7t  CB.1.16) 

i.e.,  Target  at  infinite  depth,  a  condition  redundant  for  the  problem 
defined  as  the  depth  is  considered  to  be  constant  at  Zq. 

vH  C^2-C22 

vii)  CJ3-C23 

=>  C^-Ca  (B.1.17) 

i.e.,  Target  at  endfire  or  far-field. 

viii)  A,2-A22 

IX )  A,3-A23 

A^=A2k  (B.1.18) 

i.e..  Target  motion  along  perpendicular  bisector  plane  of  the  line  joining 
the  two  sensors  or  target  at  endfire  or  target  at  far-field.  Finally, 
with  some  mathematical  manipulations  one  can  derive  the  following: 


I 


-8(b|2C|^-b,^C^2^  t3l3^15“®15’^13^“^^^14‘'13“^13‘-14^^®l-1^12 
-6(b^2t^5-b^5C^2)(a^^b,3-a,3b^^)♦2(a^2b,3-a^^bJ2)^b„c,5-b,5C,^) } 
Similarly, 

1D,(.I  -  IDacll  I  I^ac22"  ^ac21^ac11  '^acl2l  ^acll  *  exists. 

Also, 

ID„„I  .  -,,m(^^)5in(^53a)co>(^^ii3,3^„(^3j^2J, 

.4co>C^!!!a^)3ln(^i!2^)sln(^^23)sin(5!l^) 

The  above  two  terms  in  the  R.H.S.  are  identically  zero  iff 

i )  C|2*C23 
11 )  C,3-C23 

=>  Clk“^2k 

i.e..  Target  located  at  entire  or  far-field. 

hi)  A,3-A23 


(B.1.19) 


(B.1.20) 


B-6 


•v)  A, 2-^22 


=>  ^Xk~^2k. 


(B.1.21) 


i.e.  ,  Target  motion  along  perpendicular  bisector  plane  of  the  line  joining 


the  two  sensors,  or  target  located  at  endfire  or  far-field. 


v)  A,3-A23-n 


Vi  )  ^12^^22*^ 


(B.1.225 


i.e.,  target  at  infinite  depth,  a  condition  ruled  out  for  this  problem. 


As  before  one  can  derive  the  following  expression: 


ID^gl-  12(aj2C,3-a^3C,2)(a,^b^5-aj5b^^j-4(a,^c^3-a,3C,^){a,2b,5-a,5b,23 


-8Ca^2C^^-a^^c^2Hal3bl5“at5^l3)“3(a^5C^3-a^3C^5)(a|^b^2~al2^M^ 


-  6  ( a]5Cf 2~^12^1S ^  f  ^  2  ( a|2bt3  ®13^1^  ^  ^ ^  } 


(B.1.23) 


M 


Similarly, 


•  I 


^ab22  ”  ^ab2l  *^abll  ^abl2  ^  ^  ^abll  exists 


where , 


.4co,*,co5(^U^a)„„cE«;£a),.n(^^^2J,M„C^U^) 


(B.1.24) 


The  last  two  terms  are  identically  zero  if. 

i)  ^^2*^22^^ 

ii)  ^13*^23*^ 

^  C^jj*C2|j-n  (B.1.25) 

i.e.  ,  Target  motion  along  perpendicular  bisector  plane  of  the  line  Joining 

the  two  sensors, 
ill)  ^\2~^22 
iv)  C|3»C23 

Cyj'Cat  (B.1.26) 

i.e.,  Target  located  at  endfire  or  far-field. 

V)  CB.1.27) 

i.e.,  target  at  infinite  depth,  a  condition  ruled  out  for  this  problem. 

As  before  one  can  show  that 

1 0,1,1  -  12(aj2bj3-a^3bj2Ka^^b5-a|5b^.4]-4(a^^bj3-a|3b,^)Ca|2b^5-a,5b,2) 

-8(aj2b|^-a^^b^2)(®t3b|5“at5^i33~3(a|5b|3-a^3b|5)(a^^b|2”ai2*^M^ 

-6(a,2b|5-a^5bi2)(ai^b|3-a^3b|^J  ♦2(a^2^i3“®i3^i2J^®M*^i5“®i5^M^ }  (B.1.28) 

One  can  now  combine  conditions  (1)  through  (6)  and  the  respective 


equations  to  determine  the  determinant  of 


S  uf  f  ic ient  C  ond  it  ions  of  Unident  if i ability 

From  the  above  analysis  it  is  clear  that  for  determinant  of  the  main 
matrix  to  be  zero  the  determinant  of  all  2x2  submatrices  must  be 

zero.  The  details  of  the  workout  can  be  obtained  in  [29).  Also  on 
inspection  all  the  above  conditions  reduce  to  the  following  three 
important  conditions  as  follows; 

i) 

=>  (B.1.29) 

i.e.,  Target  motion  along  perpendicular  bisector  plane  of  the  line  Joining 
the  two  sensors. 

ii) 

^  (B.1.30) 

i.e.,  Target  located  at  endfire  or  far- field. 

iii)  (B.1.31) 
i.e.,  Target  located  at  endfire  or  far-field,  or  target  motion  along  the 
perpendicular  bisector  plane  of  the  line  joining  the  two  sensors. 

B,2  Step  II-TDOA  and  depression  angle  measurements 

In  this  case  one  has  to  show  that  the  1D|  exists  where; 

A.  At  At  A^  Ac 

B,  Bg  bJ  Bs 

ID  I  *  Ct  Ct  C5  (B.2.1) 

A. t.  AtIt  AtXt  A4t4  Aetc 

B, t,  B2t2  B3t3  B^t^  Bjts 

From  Ref.  114)-[16]  it  is  sufficient  to  show  that  ID^I  ♦  ID2I  ♦  ID3I 
is  non-zero,  where  ,  O2  and  D3  are  as  defined  in  Section  II. 

Existence  of  the  first  term  is  already  shown  in  Step  I  above.  If  one  works 
through  the  evaluation  of  I  Dj  I  and  I  D3  I  one  can  show  that: 


=  A^{  1 2H(K'  -X'+T' )  -  4S  (L'-Y'  ♦?' )  -  8Q  (N*  -Z'-R' ) 

+3R(M'-V'■^Q')-6P(0'-V'♦S')♦2T(F'-G'•^H'))  (B.3.2) 

where  the  constants  can  be  obtained  from  [29).  On  closer  examination  of 
Eqn.  (B.3.2)  it  is  clear  that  the  terms  on  the  R.H.S.  are  identically  zero 
if  the  submatrices  of  determinants  ,  O2  and  D3  are  zero  which  as 
in  step  one  reduces  to  the  three  main  conditions  as  in  the  Equations 
(B.29)  (B.30)  and  (B.31). 

B  .  3  I  de  n  t  i  f  i  a  b  i  1  i  t  y  condition  for  the  three-dimensional  case  when 
only  frequency  measurements  are  considered 

When  the  FIM  is  evaluated  by  considering  the  s i ng 1 e - s e ns o r 
target  frequency  measurements  obtained  with  two  sensors,  the  target  is 
generally  identifiable.  The  sufficient  conditions  of  uni  de  ntif  i  ab  i  1  i  ty 
are  developed  as  follows:  In  this  case  one  has  to  show  that  the 
determinant  of  the  matrix  0^ .  i-4,5  is  non-zero,  where  the  I  Dj  I  is 
given  by 


ID,|  = 

an 

bii 

ai2 

*^i2 

ai3 

ai4 

bi-i 

ais 

biS 

<^11 

^iZ 

d,2 

^12 

®13 

*15 

i=4,5 

(B.3.1) 

where  the  quantities  that  make  up  the  determinant  are  obtained  from 
Appendix  A.  As  before,  one  can  develop  the  determinant  along  the  first 
column  giving  the  sufficient  conditions  for  the  I  Dj  I  to  be  zero.  The 


basic  conditions  are 


B-IO 


1) 

or 

its 

minor-0 

and. 

(B.3.2) 

2) 

b„-0 

or 

its 

minor-0 

and. 

(B.3.3) 

3) 

c„^-0 

or 

its 

minor-0 

and. 

(B.3.4) 

4) 

d,^-0 

or 

its 

minor-0 

and. 

(B.3.5) 

5) 

®lk“0 

or 

its 

minor-0 

and. 

(B.3.6) 

where  i-4,5  and  (  since  the  minimum  number  of  measurements  are 

considered.  ) 

When  the  target  is  at  endfire  and  moving  along  the  line  joining  the  two 
sensors,  then  is  some  value  say  10  kn  and  Vy=0,  one  can  show  from 

the  Eq.  (B.3.2j  that 

l-(cosBj^cosAj^)^  CB.3.7) 

which  is  true,  since  Aj|^  and  Bj^^  are  both  zero. 

Similarly,  one  can  show  that  the  Eq.  (B.3.3)  reduces  to 

cosB  jijCOsA  ji^si  nB  j^cosA  0  CB.3. 8  ) 

which  is  true,  since  Aj^  and  Bj,^  are  both  zero. 

Similarly,  under  the  above  conditions  the  Eq.  (B.3.4)  reduces  to 
tanAj,jCosBj,^cosA„j-0  (B.3.9) 

which  is  true,  since  Ajj^  and  Bj^  are  both  zero.  The  minors  of 

d.,^  and  e.,^  eventually  reduce  to  zero  as  a  consequence  of  the 
Equations  (B.3.7),  (B.3.8)  and  (B.3.9). 

When  the  target  is  at  far-field  it  can  be  considered  moving 

perpendicular  to  the  line  Joining  the  two  sensors  which  means  v,>>0  and 

Vy  is  some  value  say  10  kn,  then  the  Eq.  (B.3.2)  reduces  to 

sinBj^cosAjijCosBjijCOsAji^-O  (B.3.10) 

which  is  true  since  B,^-90®  and  Aj^-O. 

Similarly,  Eq.  (B.3.3)  reduces  to 


(B.3.11) 


•  ;  -  .  •  «  ' 


(sinB,^C03A„^)^*l 

which  is  true  since  B„^-90°  and  A„^-0. 

Also  Eq.  (B.3.4)  reduces  to 
tanA„^sinB,^cosA„^»0  (B.3.12) 

which  is  also  true  since  Aj,^  »0,  as  before  the  minors  of  dj,^ 
and  ejj^  eventually  reduce  to  zero  as  a  consequence  of  the 
Equations  (B.3.10),  (B.3.11),  and  (B.3.12).  Thus  target  unidentifiable 

if  at  endfire  or  far-field. 

When  the  target  is  moving  along  the  perpendicular  bisector  plane 
of  the  line  joining  the  two  sensors  it  is  convenient  to  consider  the 
convenient  to  consider  the  complete  FIM  and  then  evaluate  the 
determinant  of  the  combined  FIM.  In  this  case  the  target  still  moves 
perpendicular  to  the  line  joining  the  two  sensors  so  that  v^-O  and 
Vy  is  some  value  say  10  kn,  in  which  case  a,^^  &  85,^  0  and 

b^  &  bj^  ^0,  since  +  -  1  8  0“  ,  and 

thus  for  motion  along  the  perpendicular  bisector  plane  of  the 
line  joining  the  two  sensors  the  target  is  identifiable. 
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